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ABSTRACT 


The objective of this study is to investigate Vertical 
Steam Condenser design optimization, in an effort to develop 
condensers with minimum weight and volume. Candidate vertical 
condenser tubes considered were of two types: (1) external 
axial flutes with smooth internal tubes, and (2) doubly en- 
hanced tubes that include internal (coolant side) heat transfer 
augmentation. A computer program (VERTCON-2) was developed to 
be used as a preliminary design tool for selection of the 
optimal condenser that meets the specified geometrical, and 
operating parameter limits. 


The computer model considers condensation in the 
presence of non-condensable gases. Also, a method was developed 
to evaluate thermal resistance of the condenser tube wails due 


to the axial flutes. 


Thesis Supervisor: Warren M. Rohsenow 
Title: Professor of Mechanical Engineering 
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L. INTRODUCTION 


A design optimization procedure for vertical steam 
condensers is established here to develop condensers with 
minimum weight and volume. The performance comparison of 
enhanced steam condensers conducted in Reference-3 suggests 
the following methodology that provides the basis for this 
optimization scheme. 

The optimization procedure is outlined as follows: 

1. Design for the following given condenser operating 
conditions: 

a. total steam condensed 

b. condenser operating pressure 
ec. coolant inlet temperature 

d. coolant flow rate 

2. This design is subject to the following specified 

parameter limits: 
ae maximum condenser length (ene 
db. maximum pumping power Ose 
ec. coolant velocity range from a minimum velocity 

) to a maximum velocity Wires) 

3. For a given heat transfer surface minimum weight 
and volume will always be achieved for a maximum coolant flow 


) 


Sain 


velocity that results in: 
a. condenser length is equal to Le ase for pumping 


power less than P or 


max’ 
bd. pumping power equal to ee for condenser length 
less than Lax’ 
4. Thus, the optimization routine proceeds as follows: 
a. select value for coolant velocity (V) 
be. calculate overall heat transfer coefficient (U) 


using the equations from Chapter II 
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c. calculate tube length and the rest of the 
condenser geometry to finally determine total condenser 
length (L.) 

d. check if pumping power (P) is less than ee 
and L.= ay 

e. if the above conditions for P and L, are not 
met,then either change V or enhancement of heat transfer 
surface 

5. a. if the calculated condenser length LSLiax 
and oe ee change to less enhanced surface until Loe are 


ae 
or Ps a 
: } > a } 
d or L. Ee nd F<P increase the surface 
; 7 — Dw we 
enhancement until either Le L. or Ps eae 


GC, alee fe se and cag aatas then decrease V until 


either L = L,, and PsP,,, or P= Pi, and LsL,., 


Further details of this optimization procedure can be 
found in Chapter V-Computer Modeling. 

Vertical condenser tubes can provide great flexibility 
in this design optimization procedure due to the different 
possible combinations of internal and external enhancement 
configurations. Thus it 1s possible to "fine tune" the 
overall heat transfer coefficient to obtain the proper con- 
denser length. 

The procedure that has been outlined here is included 
in a preliminary design program VERTCON=-2 along with the pro- 
posed vertical condenser arrangement shown in Figure-l. 

Steam enters the condenser shell along a major portion 
of the tube length. Steam lanes distribute the steam around 
the tube bundle with radial steam inflow into the tube 
bundle towards the air cooler section and air removal duct. 
Steam condenses on the vertical tubes, and the condensate 
collects in the valley of the tube fluted surface flowing 


down the tube length. The tube support plates have a dual 
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FIGURE 1 


VERTICAL ENHANCED TUBE CONDENSER 
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purpose by also serving as condensate drain-off plates. 

The condensate is stripped off the tubes and collected on 
the tube support plates, and then it is removed to the hotwell 
section via downcomer drainage tubes. Improved condensate 
control is afforded by this condensate drainage and removal 
scheme. The tube support/condensate drain-off plates 

divide the condenser up into sections, where the full depth 
of tubes in each section see relatively fresh steam. The 
tubes will see fresh steam through the depth of the tube 
bundle since the tubes are not subjected to condensate 
indundation effects because of effective condensate removal 
and vertical orientation of the tubes. A continuous duct 
that runs along the total tube length is utilized for air 
removal. Figure-1 shows a double flow coolant water circuit 
with the inlet/outlet header located in the bottom of the 
condenser. The condenser hotwell is located below and also 
encases the inlet/outlet header. This hotwell location adds 
to the total height of the condenser. 
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IT. HEAT TRANSFER CALCULATIONS 


Nomenclature 

A area (£t* ) 

a amplitude of the flute (ft) 

Cp specific heat (BTU/1bm °F) 

JO) Fe diameter (ft) 

2 helical ridge height (ft) 

f friction factor [ a( sP/e) (D/L) (2e,/V*) 

G mass flux (lbm/ft*hr) 

g. gravitational constant (lbm nee) / dl |ae hr“) 

h heat transfer coefficient (BTU/hr ft* °F) 

heg latent heat of vaporization (BTU/lbm) 

K thermal conductivity (BTU/nr ft* °F) 

ad length (ft) 

i lead of ridge (axial distance per 360° turn) (ft) 
m operand for friction factor equation 

Nu Nusselt number 

P pressure (lof/ft*) or (in-hg-abs) 

p pigen Of Kwidgingor flute (ft) 

Pr Prandtl number 

Q heat flow (BTU/hr) 

r operand for friction factor equation 

Re Reynolds number 

RE» flooding Reynolds number for flute [= MW o/wX, | ** 
Re thermal resistance of axial flute (nrft® *2/sTU) 
Re thermal resistance of scale (hrft* °%/BTU) 

a, thermal resistance of tube wall (nrft © °s/gry) 
St Stanton number 
T temperature (°F) 

ss Tsat7twl *) 

AT. mean overall temperature difference ( ‘F) 

U overall heat transfer coefficient (BTU/hr ft* °F) 
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Dimensionless velocity | at/u* | * 
friction velocity [=y7)2,/? # 
velocity (ft/hr) or (ft/sec) 


flooding axial mass flow of condensate per flute 
(lbm/hr ) 


half-perimeter length of flute (ft) 


Subscripts 


nos & Ho oe 


sat 


fluid at the bulk temperature (°F) 
condensate 

fluid or flooding 

inside, or inlet 

nominal 

outside, or outlet 

increment, or section 

saturation 

wall 


Superscripts 


s. 


dimensionless parameter 


Greek Symbols 


a 


» 


45 a 0 F 


sg 


tr 


5 


height of the condensate in the center of the 
Biuce: Cle ) 


fragment of U,+[2 -2.5in(2e/d,) + 3.75] * 


: ZB 
dimensionless group = “5 zg. 20)" oie 
dynamic viscosity (lbm/hr ft) L 


density (lom/ft?) 

surface tension (lbf/ft) 
non-dimensional group + 
apparent wall shear stress 


Reference (6) 
Reference (4) 
Reference (1) 
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TT.A ANALYSIS SUMMARY 


The rate of heat flow Q over an entire heat exchanger 
is related to the mean overall temperature difference ATS mn 
and the total heat transfer area A by the overall heat 
transfer coefficient U. 


Q= UA 4T,, (II-1) 


From the analysis of externally ridged tubes it proves 
convenient to base the heat transfer coefficient U on the 
surface area of a smooth tube having an outside diameter 
D, equal to the diameter measured over the mid-height of 
the external flutes (see Figure-2 for external flute 


geometry). 
Un= es ee atin (II-2a.) 
Opal 
where: Di= Dy + 22 (II-2b. ) 


The overall U depends on the resistances in series; 
between cooling water and the tube wall, within the tube 
wall, and between the tube wall and the working fluid. 
Fouling resistance on either side of the tube wall can 
be combined in one term Ras 

Tube wall resistance is considered to be the summation 
of an equivalent smooth tube resistance and flute resistance. 
Tube wall resistance for a smooth (non-enhanced) tube can 
be expressed as: 

Rw? a 2n(Po/Pi Pret (11-3) 
Where Def is the reference diameter on which the overall 
U is based on. For a fluted tube the overall U is based 
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on the nominal diameter D,° As will be discussed in the 
next section , the flute resistance can be written: 


-.0242 
7 [7.985 ) + 7688(= }P 


ales 


Therefore, the total tube wall resistance for the fluted 


tube can be written: 
a 
Ry= in [Py /0i. (7.985 (a)*7688(5 ) +1) (II-5) 
Zk 


Therefore, the concept of resistances in series yields: 
pe Vie es (II-6) 
U, iS W 
he no 


The overall heat transfer coefficient UF in Equation 
(II~6) does not consider non-condensable gas effects. 
Condensing in the presence of non-condensable gases will 
be taken up in Chapter Vy. 
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FIGURE 2 
Cross Section of Vertical Fluted Condenser Surface 
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FIGURE 3a. 
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Heat Plow Restriction Model FIGURE 3c. 


Application of Heat Flow Restriction 
Model to Axial Fluted Tubes 
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IIT.8 TUBE WALL RESISTANCE 


One method to evaluate tube wall resistance for fluted 
tubes would be to use Equation-3, where Do and Oe equal 
the nominal diameter D° 


R= 1n(Dy/D; )D, (II-7) 
2K 

This formulation is not sensitive to flute geometry 
such as changes in amplitude-to-pitch (a/p)ratio. In com- 
paring Equation-7 with finite element formulation differences 
of up to 30 percent in flute resistance are observed. fThis 
could amount to 10 percent difference when considering 
total tube wall resistance for tubes made with low conduct- 
ivity materials such as titanium with thin tube wall dim- 
ensions. 

What is proposed here is to model flute and condensate 
layer as a heat flow restriction, by describing two-dimensional 
heat flow through an elliptical boundary. (Figure- 3%) 

This serves as an accurate model for a heat flow restriction 
as shown in Figure-3b. 

The procedure is to apply this same model to describe 
the heat flow through the axial flutes (Figure-3c). 

Some aspects of curvilinear coordinate systems (7) 
will be reviewed in order to develop nomenclature to be 
used in later equations. 

Rectangular coordinates x,y,Z, can be expressed in 
terms of new coordinates Up» Ups Us by the equations: 


(see Figure-4) 
AS X(Uy > Uy» Us) 
y= y(u,, U5» uz) 


2= Z(u, >» Uy » us) 





a= 


The coordinates are chosen such that Uz» Us: Us» are 
orthogonal curvilinear coordinates. From examination of 
Figure-4+ it can be seen that if a particle moves from 
point P in such a way that u, and u, are held constant 
and only uy varies, then a curve in space is generated. 


Let r represent the position vector of a point P 
in space: 
= is “A “ 
r= xl ¢ yj # zk 
a 4 A : : 
where 1, Jj, and k are unit vectors in the x, y and z 


directions respectively. 
Then a tangent vector to u, curve at point P is given 


by: F 
= - dS 
a = (II-8) 
ou, 98, duUy 
Sy," are length along the u, curve 
dr is a unit vector 
0Sy 
thus: U,= nyu, 


where u,is the unit vector tangent to u, in the direction 
of increasing arc length, 


thus : h,= qSy which is the length of U, 
du 
a 


Similarly for the other coordinates: 


ae li=9) 
hi= a ae es | dS; = h;du; BE 1,2 7) 


ak 





Soe 





Sie 
ee 


a oe 
© @ e 


FIGURE 4 
Curvilinear Coordinate System 
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The unit vectors u,, Uy and U, form a right-handed 
orthogonal curvillinear coordinate system. The vectors 
Ujdu,, Up du,, and U,du, are mutually perpendicular vectors 
having as their length the arce-length differentials dsj» 
dS5» and dS... Therefore, an element of volume can be 
determined from the following vector product: 


(II-10) 
qdV= hy hoh,du,dupdu., 
Also, a vector element of surface area dA, on the surface 
u,= constant is given by the vector product: 


dA, = hohodu,du. (II-lla) 


and similarly: 
dA5= hgh, du,du, (bd) 


dAs= h hodu, du, (c) 


GENERAL PROBLEM 


Consider steady one-dimensional heat transfer through 
a thick-walled cylindrical shell a section of which is 
displayed in Figure-5. 
First assume rate of heat transfer is constant at 
every section. 
= constant 
Qa KAGT (II-12) 
ds 


Where S is defined in Figure-5, integrating Equation (II-12) 
between inner and outer surfaces with the assumption K is 
constant yields: 


T_ =T 
Le =e (II-13) 


Q= 
(1/K) 24 R 
s/A(s) 


where R is heat transfer resistance. 
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thus: R= 4f = 1 (_ds_ (II-14) 

Q KJ dA(s) 
Equation (II-14) is only applicable when heat flux is con- 
stant across the heat transfer area. In many geometries 
the distribution of heat flux across the heat transfer area 
may not be uniform, but equation (II-14) can be expanded 
for variable heat flux (16). 

Equation (II-14) can be rewritten: 


aT ds.ds (II-15) 
= -K/ sy ds, e°3 
where: dsj= h,du, 


SS meee 
ds, h,du, 


thus: Ata Sac? 
Sj hy du, 


and substituting into Equation (II-15): 


__, at If hong du,du 
f= wh du,J/, ~ h 3 


1 
thus: (u ) 
Rs ate 2 f Ltt au (rr-26) 
Q «K ff dupdu, 
ny 
(uy), 


Equation (II-16) is applicable to variable heat flux dis- 
tribution over the area of heat transfer. But this equation 
is restricted to one-dimensional, steady state analysis, 
where temperature is a function of uy only. 

Equation (II-16) can be applied to the model of two 
dimensional heat flow through an elliptical boundary. 

Referring to Figure-3 and Equation (II-16) the following 
expressions are obtained: 





A Gyo)e 


7;= constant temperature lines 


4= 
7,¥,2 coordinates, where z coordinate is 
into page 
h;= or | where Y = xi ¢ y} zk 
Sule 
i 


thus : h,” = (Sa; ”] +a) ee) (II-17) 
ou; 


in this case: 
X= c cosh” cosy 
y= ¢c sinh” siny 


from Equation (II-17): 


he= al 


substituting into Equation (II-16) 


R= 


—— of ii Er 


and: R= 3 (II-18) 
Kr 
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Merce Osn mcos / 
at ¥= 0 

X2 bde2c coshy 
and similarly: 
y= c Sinhy siny 
at ys 7/2 

Y= a= c sinh» 


therefore: 


7,2 sinh (2) 


and from a fundamental identity: 
2 
a +a +1 (II-19) 
1,5 in|? C | 


For axial fluted tubes Equations (II-18) and (II-19) 
will take the form: 


R= _? (II-20a. ) 


r= anfe@yeraye] 


where C is defined as the flute resistance coefficient. 

Equation (II-20) is in good agreement with actual physical 
situation: 

If aso this implies ?= 1n(1)s0, thus Reso 

1. Therefore a flute of zero amplitude has zero thermal 
resistance. 

2. As flute amplitude increases the flute resistance 
also increases. 

Equation (II-20) can be simplified as follows: 

1. Amplitude-to-pitch (a/p) ratios for fluted tubes 
are in the range of .1l-to-.5. 

2. The flute resistance coefficient (C) will have an 


approximate value of C# .5 thus, 





oye 


é = 
© (2) +1 =1 


and equations (II-20) simplifies to: 


R= in [c()-+1| ree 
eam 


To express thermal resistance in units of (nrft~ F/BTU ) 
multiply R from Equation (II-21) by a reference area (A Mie 
This reference area should be the same reference area 
that the overall heat transfer coefficient is based on. 
thus; 


ref 


Ra Aree (4T) (pret? R/gru) (11-22) 
Q 
a 
Roz in[c(#)+ 1]0,.5 (II-23) 
f 2K 


Values for flute resistance are calculated from a finite 
element formulation using the ADINAT (2) finite element pro- 
gram. The flute resistances obtained from ADINAT are sub- 
stituted into Equation (II-23), and the values for flute 
resistance coefficient (C) are determined at different 
values of a/p, and for various tube dimensions. One result 
is that Cc is a function of both flute amplitude a, and 
a/p ratio. Thus C as a function of a and a/p was correlated 
using multiple non-linear regression (15) analysis. 

Then a final formulation for flute thermal resistance 
was established in the following form: 


by 2 
Rez 12 I, (a) ¢ +1 Pret 
2K (II-24) 
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where C= 0, (ay (2)? 


Do» b,, b, are correlation coefficients resulting from the 
multiple non-linear regression. 

The convection heat transfer coefficient along the 
perimeter of the flute is required for input into ADINAT. 
The distribution of convection heat transfer coefficient 
was estimated by first calculating the condensate film 
thickness (5) along the perimeter using a procedure estab- 
lished by BARNES (4). The convection heat transfer at a 
particular point along the flute perimeter can be approx- 


imated by: 


h < (BTU/nrft* °F) 


c= 
where K is the thermal conductivity of the condensate, and 
is the condensate film thickness. 

As stated before, the finite element results, along 
with regression analysis were used to determine the value 
of Cc shown below: 


b dD 
Cz bo(a) (2) * (11-25) 


bow 7-985 
b,= eerste 
Do =-= 1 . 0242 


Finally, equation (II-24) can be written: 


R= in[7.985(a)*7089(8)- °°" 43 I (II-26) 
f= ig 


The final results of this analysis is depicted in 
Figure (6). Curve-1 represents flute thermal resistance 
as a function of a/p ratio for constant amplitude flutes. 
In Curve-2 both flute amplitude and a/p ratio are varying, 
but flute pitch is held constant. 
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Also, compared on this plot are the flute thermal resistance 
calculated from the finite element results, from Equation 
(II-7) based on nominal tube diameter, from Equation (II-26) 
based on heat flow restriction model. For ease of comparison 
thermal resistance have been normalized on the maximum 

range of thermal resistance considered, which is denoted 
Rp*. Figure-16 shows that the discrepancy in value between 
the finite element solution and result based on nominal 

tube diameter from Equation (II-7) increases with increasing 
flute amplitude and a/p ratio. Both Curve-l and Curve-2 
together cover a broad range of flute amplitude and a/p 
values to demonstrate the fit of Equation (II-26). 

From the results that were depicted in Figure-6 indicate 
that Equation (II-26) accurately models the flute thermal 
resistance. This equation is based on finite element sol- 
ution which must also be checked to ascertain whether it is 
a valid model. If the finite element model is developed 
methodically from a simple configuration where the results 
are known, to the more complicated model, then confidence 
can be placed on the results. In this case accurate spec- 
ification of convection heat transfer coefficients was 
important to the accuracy of the finite element solution. 
Also, a simple check can be made to determine if the finite 
element grid has the proper number of elements to model the 
problem. This is done by checking the heat flux approaching 
a boundary node from different elements. If the grid is 
sufficient there should be no jump in the value of total 
heat flux across the boundary at a specified location. 

These general guidelines were followed in modeling tube 
wall resistance for axially fluted tubes. 





ge 


Curve-l - Constant a (as .002 ft.) 

Curve=2 = Constant p (p= .00625 ft.) 
ee Finite Element Solution 
—--------- Equation (II-26) Heat Flow Restriction Model 
Equation (II-7) Average Flute Resistance 


R= 4,03 x 107 (nrrt*°r/Btv) 
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FIGURE 6 
Comparison of Flute Thermal Resistance Calculations 
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LL.C WATERSIDE HEAT TRANSFER COEFFICIENTS 





1. Smooth Internal Tube 

For the smooth internal tube the McAdams correlation 
was used for determination of the heat transfer coefficient 
for cooling water. 


Of 0.4 
ee E Onees cs | (II-27a.) 
Ky De k 
or: K 0.8 0.4 ne 
hp= 0.023 —® [Rep] [Pr] p ti la 


where: 
(a) all fluid properties are evaluated at the bulk 
fluid temperature; 
(b) 2300 Re, 10’ where; 
Rep= Reynolds number based upon hydraulic diameter 
(c) 0.5 Pry 120 where; 


Prys Prandtl number based on bulk temperature 


Z2. Helical Internal Ridging 


Internal tube augmentation was also investigated, with 
the application of integral multiple helix ridging (13) as 
shown in Figure-7. 

Friction factor data was correlated by means of the 


following equation: 


=> ee (T1-28) 
2.46 In| r4(7/Re)™| 


For the roughness of the helical internal ridge both r and m 


are treated as variables, and are tied in with tube geometry. 
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Both m and r vary with the dimensionless parameter e/1, 
where e is the ridge height and 1 is the lead of the ridge. 
In reference (13) WITHERS has made a distinction based 
on the criterion p/d = 0.36 in correlating the friction 
behavior to the internal geometry. It has been proposed 
that a shift in flow behavior occurs at p/d = 0.36. 
For higher values of v/d a greater degree of swirling could 
occur, compared with cascading of flow that takes place if 
p/d is less than 0.36. 

The heat transfer correlation equation developed from 
data for tubes of various configurations, and solved for 
the Stanton number becomes: 


Vz/8 
Ver [(e/d, )Re/f/8 


where: y= -(2.5 In(2e/4; ) + 3.75] econ) 


St= (II-29a.) 






5.68(e/p) > etn, Fh 


Then, he becomes: 


hpe K Re Pe G (IT-30) 
d 


at 


These equations are applicable to Reynolds number 
range 10,000-120,000 and Prandtl number range 4-10. 

See Appendix=-B for tube data, and friction factor 
characteristics of multiple-nelix internal ridged tubes. 
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II.D CONDENSATE HEAT TRANSFER COEFFICIENT 


In reference=-2 BARNES developed an equation for the 
average value of condensate heat transfer coefficient (h.) 
for an externally fluted tube of length L. In this 
formulation ha depends upon physical properties of the 
condensing fluid geometric factors, and mass flow rate of 


condensate in the flutes. 
.2307 %. -9226r,3 . LENO, 
a (Nu,0* ) E poli. a 
D pat 


( 11-31) 


7 h.W 0074 
Aa ~ 6027] fg | 
a Mya 


The non-dimensional group Nu,o* defined in Reference-2 is 
a function of flute amplitude-to-pitch (a/p) ratio. It 
should be noted in Equation- 31 the parameter L is the 
length of the tube between condensate drain off plates. 

The parameter W, is the flooding axial mass flow of 
condensate per flute (LBM/hr) which is also a function of 
physical properties of the fluid, and tube geometry. In 
Reference-7, PANCHEL and BELL defined the following flooding 
Reynolds number for condensate flow in the axial flutes: 

Reps 4wt ( II-32) 
wAy 


“— 


and an additional non-dimensional group: 


Z a .\4 

aes 

5 — g,(%o) ( II-33) 
- ge 

where “om 2a in the case where the flute is flooded. Also, 

the following correlation was developed for A, based on the 





= 


half-perimeter length (X,) of the flute: 
A nm 36(a/p)exp(3.33 a/p) Re, (II- 34) 


Substituting equations (32)and 33)into equation (34) and 
solving for We yields: 


Wee 8 e0” arp exp [-3-33(a) (11-35) 
9 CT P 


where We is now based on full perimeter length of the flute. 


The non-dimensional group Nu io” can also be defined 


by the following correlation: 
af 


Nu,o* = bot Sb, (a/p)* (II-36) 


where b.» b, are the correlation coefficients: 
bo= 3g (etedl 2K! 
by = 6.4526 
boa 15.265 
d35 16 .14543 
by =~6. 561166 


Equation (II-36) was obtained from multiple linear 
regression, using techniques of Reference-15. 
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PiatemorosOLbON Om PX TeRNAL ENRANCENENT CECMETRY 





The heat transfer performance of various vertical tubes 
with different levels of external enhancement were examined. 
This was done to provide guidance in the selection of external 
enhancement geometry. The basic parameters that need to 
be specified are: 

1. Flute amplitude (a) 

2. Flute amplitude-to-pitch (a/p) ratio, 
and thus flute pitch (Pp) is also set. 

3. The ratio of condensation rate-to-flooding 
condensation rate (W/W p) » 

The following plots of da ye versus a/p ratio (Figure-8) 
illustrate the relationship between these basic parameters. 
The variable al is the condensing heat transfer coefficient, 
which is normalized in the plots by the maximum value of 
lay (denoted a The curves of Figure-8 denote lines of 
constant flute pitch (p) where curve-l1 is for W/N a= Or 
and curve-2 is for WAT a= 0.25. It can be observed that 
greater enhancement is achieved with larger values of a/p 
ratios, and smaller values of flute pitch. In Figure-9 
a line of constant flute amplitude is also plotted at W/W pal. 
Figure-9 shows that increased enhancement is more readily 
achieved from increasing a/p ratio by reducing flute pitch. 

Increased steamside enhancement cannot be achieved 
Without paying a cost. The cost in this case is increased 
tube wall resistance. In Figure-10 normalized flute resistance 
for constant amplitude flutes (curve-l, az .0125") and 
constant pitch flutes (curve-2, p= .025") are plotted along 
with heat transfer data for i ol. Pigure-1ll is a 
Similar plot except curve-l is for az .0l1, and curve-2 is 


for p= .050. Flute thermal resistance was calculated using 
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2 Constant a (h _/h *) 


2 Constant p (n= /n*) 
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Equation (II-26). From Figures 10 and 11, the following 
conclusions can be inferred: 

l. When comparing on the basis of flute pitch 
at the same a/p ratio, it can be seen that flute resistance 
is larger and heat transfer enhancement is less for 
increased values of flute pitch. 

2. When comparing on the basis of flute amplitude 
at the same value of p, it can be seen that flute resistance 
is larger and heat transfer enhancement is greater for increased 
values of flute amplitude. 

This indicates that the best tube geometry is small 
amplitude flutes with larger a/p ratios, and thus smaller 
values of flute pitch in comparison with other tubes. Thus 
the benefits of increased heat transfer enhancement are ob- 
tained with minimum increase in flute thermal resistance. 

As demonstrated by Figure-8, greater heat transfer 
enhancement is achieved at smaller values of WA, but there 
is a lower limit in practical condenser design for W/W... 
Reduced values of condensate drainage plate spacing are 
required to achieve smaller values of W/W,. Close drainage 
plate spacing will complicate condenser fabrication and 
maintenance. Therefore acceptable values of condensate 
drainage plate spacing should be used as criteria for specify- 
ing W/W, in a practical condenser design. This criteria is 
used in the program VERTCON<-2 to specify ee This effect 
of Wht. on condensate drainage plate spacing can also be 
observed in Figure-8, where the drainage plate spacing at 
meant A (WAN ox .25) is equal to 1.632 ft. while at Point B 
(WAN y= Pemtneoespacine 150.622 ft. 
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TII1. CONDENSATION IN THE PRESENCE OS NON-CONDENSABLE GASES 


eg CS a Ta eT EO ——————eeue_u_eee 


Nomenclature 

A Flow area (m), (£t*) 

a Constant defined in Equation (III-6) 

B Term defined by Equation (III-l4a) 

Cor Specific heat of cooling water (J/kg°K), (BTU/lbm -F) 

Do External diameter of tube (m), (ft) 

1h Nominal diameter of axial fluted (m), (ft) tube 

D Coefficient of vapor diffusion in gas (kg/m s)/Pa, 

P (Lom/ft s)/Lot/ft2 

G Vapor mass_velocity approaching tube row (ke/m“s), 

(Lbom/ft<s ) 

g Gravitational acceleration (m/s“), (¢t/s*) 

H. Specific enthalpy of vapor (J/kg), (BTU/Lb) 

he, Specific latent heat (J/kg), (BTU/Lb) 

Specific enthalpy of vapor (Wace), (BTU/LD ) 

L Length of condenser (m), (ft) 

N Number of tubes per row 

ea Partial pressure of air in bulk of mixture (ee (1b/ft*) 

is Pressure of mixture OR (1b/ft*) 

Pe Partial pressure of vapor at condensate surface 
@a)o (lose) 

P. Partial pressurg of vapor in bulk of mixture 
Cn nytt) 

q Heat flux (W/m), (BTU/nrft*) 

Te Temperature of coolant Gee (°R, F) 

Tom Mean temperature of coolant in tube aree (°R, F) 

Ne Temperature of condensate surface (°K), (°R, °F) 

i Temperature of vapor (°K), (°R, F) 

Z Temperature of outer tube wall (°K), (°R, F 

UL, Heat transfer coefficient for convection and tube > 
resistances as defined by Equation (II-6) (BTU/hrft* °F) 

Us Overall heat transfer coefficient between vapor and 


coolant (W/m@°K), (BTU/nrft2 °F) 
V Specific volume of vapor (m?/kz), (£t?/1b) 
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We Flowrate of air (kg/s), (1bd/s) 

WT Rate of condensation per row (kg/s), (lb/sec) 
nee Flowrate of cooling fluid (kg/s), (1b/s) 

W. Flowrate of vapor (kg/s), (l1b/s) 

Ws Condenser steam flow (Kg/hr), (Lb/hr) 


Subscripts 
Subscript indicating nominal tube diameter 


Subscript indicating vapor 
Subscript indicating conditions at coolant inlet 
Subscript indicating conditions at coolant outlet 


mhrd 2 


Greek Symbols 


Cy Volume ratio of air to mixture 
=. Weight ratio of air to mixture 
Cwo Weight ratio of air to mixture at outlet of air 
cooler 
La Absolute viscosity of air (n/m s) 
am Absolute viscosity of mixture (n/m s) 
Pe Absolute viscosity of vapor (n/m s) 
é Term defined by Equation (III- ) 
v Coefficient of mass transfer (ke/m*s)/P, 


Dimensionless Groups 
Nu, Nusselt number defined by Equation (III-6) 


Re, Reynolds number of vapor gas mixture 
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HEAT AND MASS TRANSFER EQUATIONS 


A significant decrease in condensing heat transfer 
coefficient (ho) can result from the presence of very small 
amounts of non-condensable gases. The motion of the vapor 
with entrained gases moving towards a heat transfer surface, 
results in an accumulation of non-condensable gases. A 
vapor partial pressure gradient is required for the vapor 
to be able to diffuse througn the gas blanket. This results 
in a liquid-vapor interface temperature below the temperature 
of the main vapor gas mixture. The heat and mass transfer 
equations developed by CHISHOLM (6) were applied in the 
VERTCON-2 condenser design program in predicting the con- 
densation of steam in the presence of non-condensable gases. 
Diffusion coefficients obtained from experiments were used 
in this analysis by CHISHOLM (6). 

The motion of the vapor towards the condenser surface 
is described by: 

Se ACIS, (Gini) 
The rate of heat transfer per unit area is given by: 


qs Go Cn (III-2) 
and, if its assumed that only laten heat is removed, then: 


qs GN, (III-3) 


From Clapeyron's approximate equations assuming that 
a is small, then: 


Sl (III-4+) 
(TV, 





6) 


By combining equations III-l, III-2, and III-4 results 
in the governing heat transfer equation within the vapor- 


gas mixture shown below: 
a vn, 2 
Tv, (TyFs) (III-5) 
The coefficient of mass transfer is evaluated from: 


0.5 ,.-1/3. -0.6 
Nup= aren é - (III-6) 








_ ¥D 
where: Nup= -- (III-7) 
Dd 
W__,W D 
Re,= ——— 2 Geer 3) 
em 
_ FP -p 
ae —— (III-9) 
m 


For Re, 7350 


as 0.52 for first tube row 
- 0.67 for second tube row, and 
= 0.82 for third and later tube rows 
For Re < 350 
a= 0.52 for all rows 


The absolute viscosity of the mixture is determined from: 


my = [(1-E Jay # 1.6.6 4, |/(140.61€) — (III-10) 
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The coefficient of steam in air is obtained from: 


D2 2-93x10~7 fend SoS Camicalily 
Pp ia 273 


(kg/m s)/Pa 


The ratio of air volume to total mixture volume is obtained 
from: 
55 e/a (III-12) 


and the partial pressures are evaluated from the equation: 





“a 140.622 _§y 
Py” ae 


Combining equations III-5, III-6, III-7, and III-9 yields: 


2 
of 


a(t =. )?/ (III-l4a. ) 


Oa5 

D Re 2/3 

where: = p m 1/3 a bys = 

7 a Ere a ics Neg III-l14b. ) 


e) 


The heat transfer coefficient between the condensate 
and the cooling water is evaluated from: 


q= U(T,-T,) (III-15) 


where U is the overall heat transfer coefficient defined in 
equation (II-_). 
Combining equations III-14 and III-15 yields: 


3/2 
g an Ch Ty°Tl, (III-16) 


es 
2 


or: aq: (Py*e)/ tas, + a | (TII=17) 
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For given values of B, and U and temperatures T,, and 
Tr this equation can be solved for q. Equation (III-17) 
will converge for any positive initial value of q. For 
the first tube row and first iteration assume q is zero 
and the coolant temperature is the inlet value, and for 
the next row and first iteration assume gq and mean temp- 
erature of the coolant Tom from the previous row. 

In evaluating U the condensate properties are first 
evaluated at the vapor temperature. After obtaining q on 
this basis, the overall heat-transfer coefficient between 
vapor and coolant is calculated from: 


Opa) d/ (et) tii—13) 


Using the overall heat transfer coefficient and log 
mean temperature difference for the coolant, Equation 
(III-18) can be written: 


Ae eld 
Ty Tt ¢0 
substituting: Q= WpCo(T po Tp) (III-20) 


results in: 


“vf = exp(-U_7D,,LN/W.C (III-21) 

7 pT. = CXPA-UQmDyEN/W 6C, 2 

5 ee es Y's 

and the mean temperature of the coolant is evaluated from: 
Won Thm (Teo-T py WC p/U STD yLN (III-22) 


For successive iterations T is used as the cooling 


LM 
fluid temperature (T,) in Equation (III-17). At this point 





ei3 = 


the outer tube wall temperature and surface temperature of 
the condensate can be calculated, so that a more accurate 
value of U can be determined. This procedure is repeated 
with more accurate values of U, B, and Tom until successive 
values of q are within specified tolerances. 

Then, the total condensation for a tube row is cal- 


culated from: 
Ms ATDyLN/N pg (III-23) 


Standard design methods use either the BEAMA (17) 
code in the UK or the HEI (18) code in the US. Both codes 
recommend that the air ejector should be sized to reduce 
the vapor gas mixture leaving the condenser to 7.5 F below 
the saturation temperature corresponding to condenser 
pressure. As an example a vacuum of 29 in.Hg corresponds 
to a ratio of non-condensable gas to mixture flowrate (oa) 
to the ejector of: 


ewe — = 0.31 

VA 

Figure- 12 can be used to determine the flowrate of non- 
condensable gases (Wy) to be used in the sizing of the air 
ejector. In this study the BEAMA code was correlated using 
multiple linear regression (15), resulting in: 


Wyz 10.0 4 .1823x1074w_-.975x1071y 2 (qr1~-2u) 
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IV. PRESSURE DROP CALCULATIONS 


Nomenclature 
Ap Minimum Plow area between outer tubes of a tube 
pank (ft) 


A Constant defined by Equation (IV-1) 
B Constant defined by Equation (IV-1) 

C Conversion constant (1b/Ft* to in.Hg) 
D 


A Equivalent diameter of flow area (the ratio of 
flow area to wetted perimeter multplied by four) 
(ft) 


D Inside tube diameter (ft) 

D Outside tube diameter (ft) 

3 Internal tube helical ridge height (ft) 
E Dimensionless relative roughness = /D; 
g Gravitational acceleration (ft/sec) 

L Perimeter of tube bundle (ft) 

Ly Condenser tube length (ft) 


Ly sho Flow length along tube bundle perimeter defined 
in Figure-13 (ft) 


La rly 

MA Equivalent number of tube rows 

p Pitch of internal tube helical ridging (ft) 
IE Tube pitch (ft) 

ge Total pressure (in.Hg) 

i Specified condenser pressure (static) (in.Hg) 
AP, Pressure loss through tube bank (in.Hg) 

APL in Mean penetration pressure loss (in.Hg) 

AP am Mean total pressure loss around tube bundle (in.Hg) 
AP on Entrance pressure loss (in.Hg) 

me 3 Pressure loss along flow length Lo » Lz, Ly 
AP), respectively (in.Hg) 

AP Waterbox inlet loss (ft.of water) 


WI 





AP 
AP 


Bess 
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Waterbox outlet loss (ft. of water) 


Tube end loss (combined inlet and outlet), 
(ft. of water) 


Flow velocity of cooling water (ft/sec) 
Flow velocity of steam (ft/sec) 


Greek Symbols 


e 
p 


Tube surface roughness (ft) 


Steam density (lb/ft?) 


Steam viscosity (lb/hr-ft) 


Dimensionless Groups 


£ 
if 
Pon 
Re 


ae 


Friction factor 

Friction factor defined by Equation (IV-19) 
Friction factor for internally enhanced tubes 
Reynolds number 
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The total tubeside pressure drop is calculated from 
the summation of the following pressure drop components: 

1. pressure loss in condenser tubes 

2. tube end pressure loss (inlet and outlet) 

3. waterbox pressure loss (inlet and outlet) 

The pressure loss for smooth internal vertical tubes 
is evaluated using a friction factor relationship developed 
by CHURCHILL (19) as shown in Equation (IV-1) which is 
applicable for any flow regime as a function of Reynolds 
number and relative roughness of the tube. 


fe 8.0[(8/Re)*“+ 1.0/(atp )3/2 ]9-9833 (qy-a) 


where: 


Az (2.457 In(1.0A7/Pe)°°40.278)| 1 
B= (37,530/Re)+© 


= surface roughness 


5x107° r+. recommended value (20) for smooth tubing 


ty 
il 


o/Ds dimensionless relative roughness where 
D; is the inside tube diameter (ft. ) 


Pressure loss for multiple-helix internal ridged tubes 
is calculated using friction factor (f ) correlations 
evaluated from Equation (II-7) resulting from manufacturer's 
test data obtained from Reference-13. The friction factor 
data was correlated with tie as a function of Reynolds 
number, and the m & r operands with the following result: 






aioe 


2 
= 8 | — (IV-2) 
. | 2.46 ln crane 





Tube properties and friction factor characteristics for 
multiple=-helix internal ridged tubes can be found in 
Appendix<B. 

Tube end and waterbox pressure losses for single and 
double pass condensers are calculated using pressure loss 
correlations evaluated from multiple non-linear regression 
using HEI standards (18). Pressure drop data from the HEI 
code can be found in Appendix=-B. The resulting correlations 
are shown here: 


Single Pass Condensers 
Waterbox Inlet Loss: 


O04 
Waterbox outlet Loss: 
\ 2558 


Tube End Loss (Combined Inlet & Outlet): 


Pip ; 0178V." (Iv-5) 


Double Pass Condensers 


Waterbox Inlet Loss: 
O4 


Pyy= .01422V°" (Iv-6) 
Waterbox Outlet Loss: 

Pugs +00635V.°"* (1v-7) 
Tube End Loss (Combined Inlet & Outlet): 

Ppt +0356V." (IV-8) 
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Vapor condensing temperature is a function of vapor 
pressure, which follows a specific saturation temperature- 
pressure relationship. With the pressure drops encountered 
in the steam flow from the steam inlet to the air cooler 
section through the tube bundle, there is a resulting re- 
duction in condensing temperature in the direction of flow. 
This temperature reduction will affect the rate of heat 
transfer in different sections of the condenser. 

In Reference-6 HARRINGTON has noted that the condenser 
performance standards of the HEI (18) are based on static 
pressure, but that the condenser tube more closely senses 
a temperature related to the total or stagnation pressure. 
The method for evaluating the steam pressure distribution 
established by HARRINGTON along with the assumption that 
condenser performance is based on stagnation pressure was 
used in this study. Also, this can be considered as a one- 
dimensional design method where the vapor flow distribution 
is represented by a Single-flow path. The flow path is 
considered as cross flow, radially inward from the tube 
bundle perimeter to the air cooler section of the condenser. 

The pressure loss incurred by the steam flow can be 
considered as two components: 

(1) steam distribution loss for steam flow from the 
inlet to around the tube bundle 

(2) steam penetration loss for steam flow through the 
tube bundle from the tube bundle perimeter to the air cooler 


section of the condenser. 
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Steam Inlet 
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FIGURE 13 


Steam Distribution Pressure Zones 
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Steam Distribution Loss 

Refer to Figure-13for the following discussion. The 
total or stagnation pressure is calculated for pressure 
zone=-l and is evaluated from: 

py 2 Pan el Ya. (IV-9) 
g 

The steam velocity required in Equation (IV-9) and for 
the entrance loss calculation is assumed to correspond to 
the recommended maximum main steam lane entrance velocities 
as shown in Table- 1. 


Table- 1 

Recommended Steam Design Velocities 
Condenser Design Recommended Maximum Main Steam 
Pressure (in.H Lane Entrance Velocity, fps 

iL 500 

Zz 400 

3 300 

Ly 250 

5 200 


The entrance flow area will have to be sized accordingly 
with respect to this entrance flow velocity. All pressure 
losses around the tube bundle perimeter are based on the 
total pressure in zonee-l. 

The next step is to calculate the entrance pressure 
loss (Ly) from pressure zone-l to the main steam inlet lane 
as follows: 


EO Nee 
wees oes (IV-10) 


2¢g 


The next step is to calculate pressure loss due to 


iercvuLon along L,- Calculate the steam velocity from the 





ee 


flow area in this zone. It assumed the net steam flow is 
reduced proportionately along the flow length Li» Los L3, 
and Ly 
Reynolds number can then be evaluated from: 
pD.V- 
pe 


Res (y=115 


where D, is the equivalent diameter of flow area, which 
is the ratio of flow area to wetted perimeter multiplied 
by four. 

Finally the pressure loss due to friction along Lo 
is evaluated using the Tanning equation as follows: 


2 
p, - ‘tCpfloV. Giry=12 ) 
LZ ni aa a eee 
D,78 


The friction factor f is obtained from Figure-14, or 
by evaluating the following correlation obtained by using 
multiple non-linear regression analysis (15), for an assumed 
roughness value. 

f= .0436 Re" 186 (Iv-13) 


This procedure is repeated to determine the frictional 
pressure loss for L., and Ly This calculation scheme can 
be conducted with a greater number of smaller length segments 
to obtain a more accurate pressure distribution around the 
tube bundle. 

The mean distribution pressure loss around the tube 
bundle can then be calculated as follows: 


Bos 1 Fen * bo Po + by Pa + by Py (IV-14) 
a re a FL, +L, 


This mean total pressure identifies the condensing conditions 
at the tube bundle perimeter. 
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Steam Penetration Loss in Tube Bank 

Calculation of Pressure loss through a tube bank is 
accomplished by replacing the actual tube bank with an 
equivalent rectangular bank that has the same perimeter 
and number of tubes. The equivalent number of rows in the 
equivalent tube bank depth is given by: 


N= ME (1¥-15) 


Next calculate the velocity through the minimum flow 
area between the tubes. The flow area is obtained from: 


Apc L 
f > (P=D, Lp (IV-16 ) 


The steam velocity entering the tube bank can be cal- 
culated using this value for flow area (A). Then Reynolds 
number 1S evaluated from: 


Res POV, 


7 (IV-17) 


The terms (P&u) in Equation (IV-17) are calculated at 
the temperature corresponding to the mean pressure at the 
tube bank entrance. 

The tube bank pressure loss equation, which is mod- 
ified for uniformly decreasing mass flow (21) 


Z 


Dz Wt" (0.5 + Ne )_C Vs (IV-18) 
5 22g 


MieanGloLion ractor et is obtained from Figure-15, or 
by evaluating the following correlation obtained by applying 
multiple non-linear regression analysis (15). 





=k 


rts 65.5Re OF (p/p, 4257 (104Re<1000) (Iv-19a.) 


2 
fle .591Re7*136 (p/p, 71-13? (1000<ReS1.0x10~) (IV-19b) 


HARRINGTON (6) has suggested that the following equation 
be used to determine the mean tube bank pressure loss: 


2 
nm vet(0. 54 Ne) Sale (IV=-20) 


Ie 
b 2 


The mean condensing pressure in the tube bank is then 
evaluated from: 
Pom= P474Pam ~4P pm (Iv-21) 


All of the heat transfer calculations are based on the 
mean condensing pressure Oot? in the tube bank. 
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Friction Factor for Steam Distribution 
Losses 
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Friction Factor for Steam Penetration 
(Tube Bank) Losses 
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VeecOnr Uren WODE LENG 


The basic procedure for heat enchanger sizing that was 
employed in the VERTCON-1 program (3) is used in the Vertical 
Condenser Sizing Program No. 2 (VERTCON-2). In addition the 
VERTCON-2 program considers condensation in the presence of 
non-condensable gases, and the program optimizes the con- 
denser design for minimum weight and volume. Other features 
of the VERTCON<-2 program are as follows: 

1. The option of smooth internal tubes or internally 
enhanced tubes. 
2. Design of single or double pass condensers. 
3. Option of submarine, or surface ship condenser 
design. 
4, Two basic configurations are available for 
double pass condensers: 
(a) Conventional return header design 
(b) "U-Tube" construction 
5. Diagnostic messages used to warn when con- 
denser drainage plate spacing is beyond specifications based 
on heat transfer performance and mechanical design. 
6. Computation of recommended values for condenser 
drainage plate spacing and maximum main steam lane entrance 
velocity when they are not specified. 


Peooram Discription 


The program VERTCON-2 is described here, and diagram- 
matically in the simplified program flow chart (Figure-16). 
The program is based on the FORTRAN-77 standard, and is 
listed in Appendix-A. Also included in Appendix-A are 
sample program input and output files. 

The program input consists of input data files, and 
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interactive input entered at the terminal. The input data 
file consists of tube geometry information for the various 
candidate tubes to be considered, and specification of 
condenser geometry, materials and operating conditions. 
The design features are selected interactively at the ter- 
minal. See Appendix-A for further details on data input. 

The program is set up for sizing just one condenser 
or a series of condensers. After inputting data the pro- 
gram prompts at the terminal for selection of design features. 
If condenser drainage plate spacing and maximum main steam 
lane entrance velocity are not specified the program cal- 
culates recommended values (Branch-A and Iteration-A, 
Figure-16). The determination of these values are based 
on heat transfer performance and mechanical design to provide 
proper tube support. If these values are specified then 
the program checks if condenser drainage plate spacing 
is beyond recommended values (Branch-B, Figure-16). If 
the values specified exceed those dictated by proper design, 
then diagnostic warning messages are printed at the terminal 
and in the program output. 

After input of interactive data at the terminal, the 
program conducts the selection of flow velocity and tube 
geometry to optimize the design as follows: 

1. For the first candidate tube geometry the 
condenser is sized at the specified minimum flow velocity 
w= V.. ie 

2.(a) If condenser length (L) is less than specified 


in 


maximum condenser length (Linax) and pumping power is less 
than specified maximum pumping power ee, then Step 1 
is repeated for a larger value of flow velocity. 


2.(b) This process is repeated until Le Bene or 


P= reas? 


Smelt eG Vin? as aes then this candidate tubde is 


rejected because of insufficient enhancement. 





oe ee 
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4, Steps 1 through 3 are repeated for all can- 
didate tubes. 

The condenser sizing in Steps 1 through 4 is evaluated 
as follows: 

1. In iteration B and C (Figure-16) the total 
effective tube length for the condenser is calculated by 
adding N tube sections of length equal to the drainage plate 
spacing until a proper heat balance is obtained. Iteration B 
adds up N-l of these sections of equal length. Iteration C 
then calculates the final section length, which may have 
a length less than the specified drainage plate spacing. 

2. After the total effective tube length is de- 
termined the condenser geometry, weight, volume, and system 
power are calculated for each condenser. 

It should be noted that the influence of non-condensable 
gases has not been considered yet, and a multiplicative 
factor (Fy) is simply used to account for non-condensable 
gas effects. 

At this point the minimum condenser sizing for each 
tube configuration are compared to select the condenser 
configuration with overall minimum weight and volume. 

If the option is selected, the program will now consider 
the effects of non-condensable gases. Using the tube bundle 
and shell geometry of the selected condenser configuration, 
the mean pressure around the tube bundle is determined. 

With the conditions identified at the tube bundle perimeter 
the program steps through the condenser tube bundle row by 
row to calculate the effective tube length at local conditions 
including non-condensable gas effects. 

If the condenser length calculated from the consideration 
of non-condensable gases differs from the original calculated 
value, then a new value for the multiplicative factor Fc 
is selected and the complete process including tube and flow 
velocity selection is repeated until these length values 





“Bio 


are equal. By conducting the calculations in this manner 
the non-condensable gas effects will be evaluated a minimum 
number of times, thus significantly reducing computer 
solution time. There is a limit to the number (Nr) of 
these iterations that will be completed by the program. 
If after Ny iterations there has not been proper convergence, 
then warning messages are displayed on the terminal screen 
and printed in the program output. This indicates that one 
of the specified condenser parameters is out of limits for 
the particular condenser design. Non-convergence would 
most likely be due to factors that influence steam-side 
heat transfer performance. One culprit could possibly be 
too small a value for tube pitch-to-diameter ratio, that 
would cause excessive steam side pressure drop. 
The tabulated program output displays for each condenser 

the following information: 

(1) All input data and design features selected. 

(2) If applicable, diagnostic warning messages. 

(3) Detailed weight breakdown of condenser 
components. 

(4) Total condenser volume, center of gravity, 
and total weight in dry and wet conditions, and system 


pumping power. 
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APPENDIX A 


DATA INPUT 
See Table-3 for a sample input data file 


INPUT DATA FILE 


1. Assigned Logical Operator - 20 to specify input 
device. See lines 47-52 of progress listing for 
further details. 


2. Free Format 


3. Variable names starting with I, J, K, L, M, or N 
are integer variables and others are real variables, 
except for specific cases as noted. 


LINE 1 
Enter: SUB or SUR 
SUB - Submarine Condenser Design 
SUR - Surface Ship Condenser Design 


LINE 2 
Enter; NUMBT ,LMAX, PMAX, VMIN, VMAX 
NUMBT=- Number of Tube configurations to be tested 
LMAX - Maximum condenser length (ft) 
PMAX -=- Maximum pumping power (hp) 
VMIN - Minimum coolant flow velocity (ft/sec) 
VMAX - Maximum coolant flow velocity (ft/sec) 
Note: LMAX is a real variable 





mers: TCL, TSAT, QSTM, STMLD, PSAT PD 
sont - Inlet Coolant Temperature ( F) 


TSAT - Steam Inlet Saturation Temperature ( F) 
QSTM = Quality of Inlet Steam to Condenser 
STMLD - Steam Condenser (lb/hr) 
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TABLE=Z3SAMPLE INPUT DATA FILE 
ee Geer hom eC ASH. 


SUB 

he? hess Ae 27S 8 es 15 

Goe elasson es. 257000.,6.5,1.35 
Dee Ghani 5625. .5625.2.64,488., 
Slates nee 488. 555.52 62.,1.5657,.56,488. 
eS G25 Ol 2245.0. ,0. 
48.,0.,0.,9.6,.00033 

PS 5 Gr Ol. .245..0204..191 
48.,.59,.00197,9.6, .00033 

ee S562, 0. 1245. .0125,.4750 
No eiG2 0 0-96. .00033 

D555) 625.001, .245..024,.0949 
48.,.58,.0075,9.6, .00033 




















on 
Si, 


LINE 4 
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- Condenser Operating Pressure (in.Hg) Absolute 


- Tube pitch-to-diameter ratio 


Emcer: TIS, 


JABS: 

TOTS 
TSHL 
OL 
THW 

THDR 
PSD 


Note: ITSD 


LINE 5 


DOS ero lcon stn,» LDR, LITSD 


Thickness of internal tube sheet (in.) 
Thickness Outer Tube Sheet (in) 
Thickness of condenser shell (in. ) 


uh 
Ab 
g 


hickness tube support plate (in.) 
hickness of hot well plate (in.) 
hickness of Headers (in.) 


Internal tube sheet density (1b/Ft?) (real variable) 


ih 


macer: OTSD, 


OTSD 
SHLD 
LoL D 
HWD 

HDRD 
TBD 

ETUBE 


TCO VER 
COVERD 


LINE 6 


s a real variable. 


Sep. loeb no, HORD, TBD, ETUBE, “TCOVER, COVERD 
Outer tube sheet density (lb/ft?) 

Shell density (lb/ft?) 

Tube support plate density (lb/ft?) 

Hotwell plate density (lb/ft?) 

Header density (lb/ft?) 

Tube bundle density (1b/ft?) 

Modulus of elasticity of tube material (lb/in“) 
Thickness of U-tube header cover (in.) 

Density of U-tube header cover (in.) 


merc: NTYPE, DL, DW, AE, AEPE, EI, PI 
NTYPE - Tube type: a. for smooth internal tube 


Dat 
DW 
AE 
AEPE 


RI 


b. for doubly enhanced tubes 
Internal tube diameter (in.) 
Tube Wall diameter (in. ) 
External flute amplitude (in.) 
Amplitude-to-pitch ratio of external ridging 


Internal helix ridge height (in.) 





eae 


ett - Pitch of Internal Helix (in.) 


Note: For smooth internal tube (NTYPE-1) enter @% for FI, 
PI. Values of EI and PI are obtained from tube 
data in Appendix<3. 


LINE 
Enter: FLUTE, MI, RI, KW, RSCALE 
FLUTE - Number of flutes per tube 


MI - Operand for friction factor, helical internal 
ridging (real variable) 

RI - Operand for Friction Factor, Helical internal 
ridging 

KW - Conductivity of tube wall (BTU/hr-ft- F) 


RSCALE - Fouling resistance (hr-ft- F/BTU) 
Note: MI is a real variable. 





To consider more tube configurations, repeat lines 6 and 7 
for each additional tube to be included in the analysis. 


INTERACTIVE DATA INPUT 
Enter data at the terminal as follows: 


When the program commences, the following message 
appears on the terminal screen: 


#% WELCOME TO PROGRAM: VERTCON, VERSION-2 ** 
ENTER PROGRAM RUN NUMBER 


Enter Run number, which can be an integer from 0-to-999. 
Next at the terminal appears: 


ENTER THE NUMBER OF CONDENSERS TO BE SIZED 
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Enter the number of condensers to be sized, which is also 
an integer from 1-to-999. 


Next at the terminal appears: 
ENTER "1" FOR SINGLE PASS, OR "2" FOR DOUBLE PASS CONDENSER 


Enter the proper integer number 1 or 2 for single or double 
pass condenser. 


Next at the terminal appears: 


TWO BASIC CONFIGURATIONS ARE AVAILABLE FOR DOUBLE 
PASS CONDENSERS : 


(1) CONVENTIONAL RETURN HEADER DESIGN 
(2) "U-TUBE" CONSTRUCTION 


ENTER 1 OR 2 FOR CONFIGURATION SELECTION 


Enter the proper integer number 1 or 2 for configuration 
selection. 


Next at the terminal appears: 


ENTER: (1) COOLANT FLOW RATE (GPM) 


Enter as real number with units of gallons per minute. 


Next_at the terminal appears: 


DO YOU WANT TO INCLUDE IN THIS ANALYSIS THE EFFECTS 


OF CONDENSATION IN THE PRESENCE OF NON-CONDENSABLE 
GASES YES OR NO 


Enter YES or Y, NO or N for proper choice. 





Behm 


Next at the terminal apvears: 


DO YOU WANT TO SPECIFY MAXIMUM MAIN STEAM LANE 
ENTRANCE VELOCITY (FT/SEC)? 


YES OR NO 


Pacer YES or Y, NO or N for proper choice. 


If yes is entered, Next at the terminal appears: 


ENTER MAXIMUM MAIN STEAM LANE ENTRANCE VELOCITY (FT/SEC) 


Enter this as a real number. 
If no is entered, Next at the terminal appears: 


PROGRAM WILL SELECT RECOMMENDED MAXIMUM MAIN 
STEAM LANE ENTRANCE VELOCITY USING CONDENSER 
OPERATING PRESSURE AS SELECTION CRITERIA. 


Next at the terminal avvears: 





DO YOU WANT TO SPECIFY CONDENSATE DRAINAGE 
PLATE SPACING (FT. )? 


YES OR NO 


Enter YES or Y, NO or N for proper choice. 
If yes is entered, Next at the terminal appears: 





ENTER CONDENSATE DRAINAGE PLATE SPACING (FT.) 


Enter this as a real number. 


If no was entered, Next at the terminal appears: 


PROGRAM WILL SELECT RECOMMENDED CONDENSATE 
DRAINAGE PLATE SPACING 





Note: At this point, entry of Interactive Data is completed 
for the first condenser. The program now completes 
the sizing for the first condenser and will prompt 
the terminal for further interactive data for any 
additional condensers. 
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TABLE 4- Properties of Multiple-Helix Internal Ridged Tubes 


py I 











Root Fin Outside 1D 
Tube Nominal diameter d,, count, area, (max) d;, No. of 
Symbol no. diameter, in in fins/in ft? /ft in starts 
A i2 3 0.623 (Stripped) 0.163 0.573 5 
> 30 3 0.627 26.4 0.640 0.575 6 
Vv 22 3 0.625 26 0.640 0.569 6 
q 27 a 0.622 26.5 0.640 0.572 6 
A 28 7 0.624 26 0.640 0.573 5 
& 31 3 0.624 26.6 0.640 0.576 3 
v. 29 3 0.625 26.1 0.640 0.575 5 
| 38 A 0.633 38.0 0.830 0.572 8 
o. 37 3 0.624 38.5 0.901 0.574 10 
40 4 0.627 2s 0.689 0.561 10 . 
© 41 3 0.628 38.1 0.852 0.572 12 
4 9 F 0.628 38.5 0.901 0.575 10 
0 21 4 0.740 26 0.640 0.684 G 
a 19 3 0.745 (Stripped) 0.195 0.692 6 
44 4 0.627 41.0 0.574 6 
= 43 3 0.628 41.0 0.573 10 
O 42 5 0.626 41.0 0.901 0.573 10 
46 3 0.627 Ae 0.573 10 
45 : 0.627 ay) 0.577 10 
Oo 13 1 0.883 26 0.841 0.820 6 
~ 32 l 0.877 26 0.871 0.825 6 
& 25 1 0.880 26 0.841 0.816 6 
e:) 24 ] 0.878 26 0.841 0.814 6 
° 23 1 0.880 26 0.841 0.816 6 
© 26 l 0.864 27.3 0.841 0.815 6 
Internal ridging 
Internal aspect ratios ID of 
Height e, Pitch p, - envelope 
Symbol in in e/d, e/p p/da; tube, in 
A 0.0125 0.475 0.0218 0.0263 0.829 1.000 
b> 0.0162 0.279 0.0282 0.0581 0.485 1.000 
Vv 0.0165 0.320 0.0288 0.0516 0.558 1.000 
re 0.017 0.385 0.0297 0.0442 0.673 1.000 
A 0.019 0.475 0.0332 _ 0.0400 0.829 1.000 
> 0.0198 0.287 0.0344 0.0690 0.498 1.000 
v 0.0207 0.469 0.0360 0.0441 0.816 1.000 
a 0.018 0.212 0.0315 0.0849 0.371 1.000 
© 0.0200 0.166 0.0348 0.1205 0.289 1.000 
0.6175 0.170 0.0312 0.1029 0.303 -1.000 
oO 0.0215 0.138 0.0376 0.1558 0.241 1.000 
<j 0.0204 0.191 0.0355 0.1070 0.332 1.000 
oO 0.017 0.391 0.0249 0.0435 0.572 1.2348 
= 0.017 0.285 0.0246 0.0596 0.412 1.000 
0.021 0.207 0.0366 0.101 0.361 _ 
0.021 0.124 0.0366 0.169 0.216 a 
O . 0.024> 0.0949 0.0419 0.253 0.166 1.000 
0.021 0.094 0.0367 0.223 0.164 = 
0.015 0.094 0.0260 0.1596 0.163 = 
oO 0.0178 6.333 0.0217 0.0535 0.406 1.5936 
® 0.0193 0.335 0.0234 0.0576 0.406 1.5936 
3 0.0205 0.340 0.0253 0.0603 0.420 1.5936 
re: 0.0205 0.330 0.0252 0.0621 0.405 1.5936 
9 0.0205 0.338 0.0251 0.0607 0.414 1.5936 
© 0.021 0.340 0.0258 0.0618 0.417 1.5936 
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TAS. 5 Friction Factor Characteristics 
of Multiple-Helix Internal Ridged Tubes 








Tube 
no. rn 4 
12 0.762 0 
30 0.64 —0.00039 
22 0.697 —0.00017 
On O72 —0.00026 
28 0.72 — 0.00028 
31 0.61 —0.00095 
29 0.68 —0.00064 
38 0.61 —0.00109 
37 0.54 —0.00392 
40 0.627 —0.00080 
4] 0.57 0.00259 
9 0.59 '—0.00197 
21 0.70 —0.00014 
19 0.626 0 
44 0.53 —0.00180 
43 0.55 +0.00017 
42 0.58 + 0.00750 
46 0.54 + 0.00275 
45 0.52 —0.00159 
13 0.63 +0.00024 
32 0.70 + 0.00090 
25 0.645 + 0.00032 
24 0.64 —0.00035 
23 0.637 + 0.00028 
26 0.64 


: +0.00018 


FIGURE 17 Friction Factor Curves for 
Multiple-Helix Internal Ridged Tubes 
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FIGURE 18a 
Operands for friction factor equation versus geometric aspect 
ratio e/l for multiple-helix internal ridging with p/d,2 0.36 
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FIGURE 18b 


Operand for friction factor equation versus geometric aspect 
ratio e/l for multiple-helix internal ridging with p/d, < 0.36 





-110- 


FIGURE =1 
WATERBOX AND TUBE END LOSSES 


SINGLE PASS CONDENSERS 
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PIGURE-20 
WATERBOX AND TUBE END LOSSES 
TWO PASS CONDENSERS 
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